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Abstract. The location of B supergiants in the Hertzsprung-Russell diagram (HRD) represents a long-standing problem in 
massive star evolution. Here we propose their nature may be revealed utilising their rotational properties, and we highlight a 
steep drop in massive star rotation rates at an effective temperature of 22 000 K. We discuss two potential explanations for it. On 
the one hand, the feature might be due to the end of the main sequence, which could potentially constrain the core overshooting 
parameter. On the other hand, the feature might be the result of enhanced mass loss at the predicted location of the bi-stability 
jump. We term this effect "bi-stability braking" and discuss its potential consequences for the evolution of massive stars. 
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1. Introduction 

The large number of B supergiants as well as their location 
in the Hertzsprung-Russell diagram represents a long-standing 
problem in massive star evolution (e.g. Fitzpatrick & Garmany 
1990). Even the most basic question of whether B supergiants 
are core hydrogen (H) burning main sequence (MS) or helium 
burning objects has yet to be answered. Here we propose their 
nature may be revealed utilising their rotational properties. 

On the MS, O-type stars are the most rapid rotators known 
(with vsini up to 400 km s _1 ), but B supergiants rotate much 
more slowly (with vsini ;S 50km s _1 ), which has been attributed 
to the expansion of the star after leaving the MS. Hunter et al. 
(2008) noted a steep drop in rotation rates at low gravities (log 
g < 3.2) and suggested the slowly rotating B supergiants to be 
post-MS. The steep drop was also used to constrain the core 
overshooting parameter a ov in massive star models (Brott et al. 
2010). The slowly rotating B supergiants are also cooler (with 
T e ff below ~22 000 K) and usin/ is observed to drop steeply 
below this T e ^. Here we introduce an alternative explanation 
for the slow rotation of B supergiants: wind-induced braking 
due to bi-stability, or bi-stability braking (BSB). 

Mass loss plays a crucial role in the evolution of massive 
stars. Whilst a large amount of attention has been directed to- 
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wards the role of stellar winds in terms of the loss of mass, as 
winds "peel off" the star's outer layers (Conti 1976), much less 
effort has been dedicated to understanding the associated loss 
of angular momentum (but see Langer 1998, Meynet & Maeder 
2003). Yet the angular momentum aspect of these winds may 
be equally relevant for understanding massive stars as the loss 
of mass itself, possibly in a mass range as low as ~1O-15M . 

We first recapture the physics of bi-stable winds and BSB 
(Sect. 12. U . before presenting the current knowledge of rota- 
tional velocities of massive stars. We note a steep drop at 
~22 000 K (Sect. [3) and propose two possible explanations 
for it. In the first one, the drop is due to the separation of MS 
objects from a second population of slow rotators (Sect. 14. 11 1, 
whilst in the second one the slow rotation is the result of BSB 
(Sect 02). 

2. The physics of the mass loss bi-stability jump 

The BS-Jump (Pauldrach & Puis 1990) is a theoretically pre- 
dicted discontinuity where wind properties change from a mod- 
est M, fast wind, to a higher M, slow wind, when the effective 
temperature drops below ~22000 K. Vink et al. (1999) pre- 
dicted an increase in the mass-loss rate by a factor of ~5-7 (see 
the blue dotted line in Fig. Q]) and a drop in the terminal wind 
velocity by a factor of ~2. Here the reason for the jump is an 
increased flux-weighted effective number of iron lines due to 
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the recombination of Fe iv to Fe m (and not necessarily related 
to the optical depth of the Lyman continuum). In fact, the tem- 
perature of the BS-Jump was found to be weakly density de- 
pendent, with the BS-Jump starting as high as -26 kK for the 
higher mass models and dropping to -22.5 kK for the lower 
mass models at ~2OM (Vink et al. 2000). 

Whilst the predicted drop in terminal wind velocity across 
the BS range has been confirmed (Crowther et al. 2006), the 
issue of a jump in mass loss is controversial. The jump may 
have been confirmed in radio data that suggest a local M max- 
imum at the predicted temperature (Benaglia et al. 2007), but 
the predicted rates on the cool side of the jump are up to an 
order of magnitude larger than those found from state-of-the- 
art NLTE models (Vink et al. 2000, Trundle & Lennon 2005, 
Crowther et al. 2006, Markova & Puis 2008, Searle et al. 2008). 
To gain more insight into this mass-loss discrepancy one way 
forward is to search for other physical effects in the bi-stability 
region, which might assist us to unravel whether B supergiant 
mass-loss rates are as high as predicted, or as low as the spec- 
tral modelling suggests. Here we outline one such approach, 
involving stellar rotation rates in the region of the BS-Jump. 
Vink (2008) pointed out that the temperature of the bi-stability 
jump coincides with the position where the rotational veloc- 
ities drop to smaller values (below 100 km/s), and suggested 
this might be due to BSB. 

2.1. Bi-stability braking (BSB) 

Employing stellar evolution models which include mass loss 
and rotation, we test whether the predicted increase in M at the 
bi-stability jump might lead to slower rotation. FigureQjshows 
the run of the mass-loss rate and predicted rotational velocity as 
a function of temperature for a 40 M star with an initial rota- 
tional velocity of 275 km s 1 . It shows a drastic drop in surface 
rotation rates for massive stars around 22 000 K, which is due 
to BSB in our models. When we do not increase the mass loss 
due to the BS-Jump, the stars remain rotating rapidly - despite 
the stellar expansion, as angular momentum is transferred from 
the core to the envelope. Bi-stability braking can only be effi- 
cient if the star spends a significant amount of time, i.e. part 
of its MS evolution, on the cool side of the BS-Jump. In our 
present standard models BSB only occurs above a critical mass 
of -30 M Q for the Galaxy, -35 M Q for the Large Magellanic 
Cloud (LMC), and -50 M for the Small Magellanic Cloud 
(SMC). In these models we employed a core overshooting pa- 
rameter a ow of 0.335 of a pressure scale-height. A higher value 
of a ow lowers the critical mass. For instance, in a test calcula- 
tion with ff ov =0.5 at LMC metallicity, BSB occurs already for 
a 20M Q star. 

3. The steep drop in rotational velocity at 22000 K 

Howarth et al. (1997) catalogued vsini values for 373 OB stars, 
with roughly half of them being supergiants (luminosity class 
I; red diamonds). We plot the vsini values of this large and 
uniformly determined data-set in Fig. [2] The figure shows a 
drop in vsini for stars hotter than 22 kK with values as high 
as -400 km/s for all objects and as high as -250 km/s for the 
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Fig. 1. Mass-loss rate (blue dotted) and rotational velocities 
for a 40 M star with an initial rotational velocity of 275 km 
s -1 , including the predicted BS-Jump (red solid) and without 



it (dashed green). A core overshooting parameter ff 
was employed in these models. 
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Fig. 2. Projected rotational velocity vsini of the Howarth et al. 
(1997) dataset of Galactic OB supergiants (red diamonds) and 
non-supergiants (blue triangles) as a function of T e ff (converted 
from spectral types using Martins et al. 2005 and Crowther 
et al. 2006). We note that vsini drops from values as high as 
~400/250km s" 1 to values below 100 km s _1 at -22 kK. 



supergiants only to values that all fall below 100 km/s for the 
cooler objects. In other words, we identify a general absence 
of rapidly rotating B supergiants^. As the stars in this data-set 
have not been analysed in detail, we resort to the results from 
the flames Survey of massive stars (Evans et al. 2008), which 
involves data from the Galaxy, the LMC, and the SMC. 

Figure [3] shows usim versus effective temperature for the 
flames data, where we again note a steep drop in vsini from 
values as high as -400km s _1 to values below 100km s _1 . The 
data selection is a non-trivial undertaking, as we wish to opti- 



1 Note that the remaining broadness in the B supergiant spectra may 
(partly) be due to macro-turbulence in addition to, or instead of, rota- 
tional broadening (Conti & Ebbets 1977, Aerts et al. 2009). 
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Fig. 3. Rotational velocities vs. for all flames objects 
with evolutionary masses above 15 M . Luminosity classes are 
shown as blue pluses (luminosity classes n-v) and red stars 
(luminosity class i). The LMC evolutionary tracks including 
the predicted BS-Jump are shown in grey with initial v rot = 
250km/s for five masses of 15, 20, 30, 40 and 60 M Q . It can 
be noted that the critical mass for BSB is ~35M Q in the LMC. 
The steepness of these tracks can be compared to the angu- 
lar momentum conservation case, drawn as grey dotted back- 
ground lines. The black dots on the tracks represent 10 5 year 
time-steps. 



mise the sample size to sample homogeneity - in the presence 
of selection effects. The reason we employ a cut-off mass at 
15M is that for the largest subset, i.e. that of the LMC, there 
is a detection limit that runs from ~2OM at the hottest temper- 
atures to ~1OM at the cool part of the HRD (see Brott et al. 
2010). For this reason, we choose an intermediate value of 15 
M as a minimum value. If we had chosen a higher mass cut- 
off, the drop feature would shift to a somewhat higher T e ff (up 
to 27 kK). If we had opted for a lower mass cut-off, the feature 
shifts to r e ff = 20 kK. Having noted this, tens of manual trials 
have shown that the sheer drop feature itself does not depend 
on a particular choice of mass range, and given the presence 
of the drop in both Figs. [2] and [3] as well as data presented in 
Fraser et al. (2010), we argue the drop feature is ubiquitous in 
the entire mass range ~1O-6OM . 

The dark grey lines overp lotted in Fig. [3] show evolution- 
ary tracks for the intermediate metallicity of the LMC with 
v m = 250km/s by Brott et al. (2010) for masses of 15, 20, 30, 
40 and 60 M . The black tick-marked dots on the tracks repre- 
sent evolutionary time-steps of 10 5 years, and are intended to 
facilitate the comparison with observations. 

4. Two possible interpretations for the drop in ( sin/ 

In Sect. [3] we highlighted the steep drop in the rotation rates 
of massive stars, but we have yet to provide an explanation for 
it. The question is whether the absence of rapidly rotating B 
supergiants is the result of BSB, or if the cooler slow rotators 
form an entirely separate population from the hotter MS stars. 
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Fig. 4. HRD of the flames survey of massive stars. See the cap- 
tion of Fig.[3]for an explanation of the symbols. 
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Fig. 5. Nitrogen abundance vs. T e g for the LMC subset. See the 
caption of Fig.[3]for an explanation of the symbols. 



4.1. The case for two populations 

The cool objects (red asterisks) in Fig.|3]and the HRD of Fig. [4] 
are supergiants of luminosity class i, whilst the fast rotators are 
dominated by dwarfs (blue pluses). Although it is by no means 
obvious that supergiants cannot be in a H burning phase, the di- 
vision in logg might imply that we have a population of rapidly 
rotating MS objects on the one hand, whilst observing a pop- 
ulation of slowly rotating evolved supergiants - which have 
somehow lost their angular momentum - on the other hand. 
Currently, we do not have sufficient information with respect 
to the evolutionary state of these cool supergiants. In princi- 
ple, this part of the HRD can be populated with the products 
of binary evolution, although this would normally not be ex- 
pected to lead to slowly rotating stars. Alternatively, one could 
envision the cooler objects to be the product of single star evo- 
lution, e.g. post-RSG or blue-loop stars, but the key point is that 
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within the context of the two population interpretation, they are 
not core H burning. 

A potential distinguishing factor between the "two popu- 
lation scenario" and BSB is that of the chemical abundances. 
We present LMC N abundances versus effective temperature 
in Fig. [5] noting that N abundances could only be derived for 
a subset of our objects shown in Figs. [3] and [4] As the LMC 
baseline [N/H] equals ~6.9, the vast majority of slow rotators 
is found to be strongly N enhanced. Although rotating models 
can in principle account for large N abundances, the fact that 
such a large number of the cooler objects is found to be N en- 
riched suggests an evolved nature for these stars. 

4.2. The case for BSB 

The second explanation for the steep drop in rotation rates is 
that both the objects cooler and hotter than 22 000 K reside on 
the MS, and that it is BSB that explains the slow rotation of the 
cooler B supergiants. The main argument for BSB is that it is 
predicted at the temperature where the rotational velocities are 
found to drop steeply. The evolutionary tracks in Fig.|4]indicate 
that the MS for the highest mass stars indeed appears rather 
broad, reaching as far as the BS-Jump temperature at 22 kK, 
and beyond. Therefore, mass loss seems capable of removing 
a considerable amount of angular momentum during the MS 
evolution for the highest mass stars. 

5. Discussion 

In principle it is possible that both effects of "two populations" 
and BSB occur simultaneously at 22000 K, with BSB occur- 
ring above a certain critical mass, and the "two population" sce- 
nario taking over in the lower mass (10-20 M ) range, but this 
situation might appear somewhat contrived. The strongest ar- 
gument for the "two population scenario" are the large N abun- 
dances of the B supergiants, whilst the strongest argument for 
BSB is that the drop is observed at the correct location (whilst 
no such coincidence would be expected for the alternative in- 
terpretation). 

Using our standard models, BSB can only operate above a 
certain critical mass and would not be able to explain the steep 
drop in rotational velocities of stars below the critical mass. 
The reason BSB does not operate at lower masses in our stan- 
dard models (of Fig. [3} is that the drop feature has been used 
to constrain the core overshooting parameter of a ov = 0.335. 
The applicability of BSB could be pushed to lower masses if 
the MS lifetime were extended. This could be achieved by in- 
creasing a ov . When we enlarge a m to 0.5, BSB also occurs at 
20M o for our Galactic and LMC models. What is clear is that 
the critical mass is model-dependent. For instance, the solar- 
metallicity models of Meynet & Maeder (2003) show BSB in 
the lower (~15-2OM ) range. 

We point out that if BSB were the correct explanation for 
the drop feature all the way down to ~1OM , we would require 
a very large core overshooting parameter, and the consequences 
would be far-reaching. For instance, it would imply that B (and 
even A) supergiants are MS objects burning H in their cores. 
This would potentially solve the long-standing problem of the 



presence of such a large number of B supergiants. Moreover, if 
BSB could work for the entire mass range, it would also have 
profound implications for the Blue to Red (B/R) supergiant ra- 
tio that has been used to constrain massive star models as a 
function of metallicity for decades. Furthermore, if the absence 
of rapidly rotating B supergiants is due to BSB, one might won- 
der what this would imply for the evolutionary state of the pre- 
sumably rapidly rotating B[e] supergiants. The rapid rotation of 
these extreme objects could possibly be related to close binary 
evolution or merging (Pasquali et al. 2000), but this requires 
future investigation. 

If BSB would indeed occur in the lower mass range (down 
to ~1OM ), one should be aware that the derived overshoot- 
ing parameter of 0.335 becomes a lower limit and that the real 
value becomes larger. Although this would be consistent with 
the suggested increase in a ov with stellar mass (Ribas et al. 
2000), such a large value of a m might be considered uncom- 
fortable, as the highest mass data-point in Ribas et al. is based 
on one binary star, V380 Cyg, for which the results have been 
challenged (Claret 2003). 

To summarise, we have presented two potential explana- 
tions for the steep drop in rotation rates at 22 000 K. Currently, 
we have insufficient information to decide which one is cor- 
rect. In any case, our study demonstrates the important role of 
mass loss for massive star evolution, and especially the impor- 
tance of specifics in its dependence on the stellar parameters. 
Furthermore, we have highlighted the significant influence of 
mass loss on the angular momentum transport in massive stars. 
Last but not least, BSB may offer a novel method of diagnos- 
ing the effects of mass loss via its influence on the angular 
momentum. Current analyses yield controversial results with 
respect to the existence of a BS jump. On the one hand, the 
predicted drop in terminal wind velocity across the BS range 
has been confirmed (Crowther et al. 2006). On the other hand, 
for temperatures below the BS-Jump, the mass-loss rates ob- 
tained from spectral modelling are generally much lower than 
predicted (Vink et al. 2000, Crowther et al. 2006). 

A simultaneous investigation of the abundances, mass loss, 
and rotational properties of a large sample of massive stars, e.g. 
with the flames n Tarantula survey (Evans et al. 2009), would 
be most helpful to settle these issues. 
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